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Abstract 

Sterile neutrinos with masses ~ 0.1 eV or higher would play an important role in astrophysics 
and cosmology. We explore possible signatures of such sterile neutrinos at long baseline experi- 
ments. We determine the neutrino conversion probabilities analytically in a 4-neutrino framework, 
including matter effects, treating the sterile mixing angles ^14,^24; ^34; the deviation of ^23 from 
maximality, as well as ^13 and the ratio AmQ/Am^^^ as small parameters for a perturbative ex- 
pansion. This gives rise to analytically tractable expressions for flavor conversion probabilities from 
which effects of these parameters can be clearly understood. We numerically calculate the signals 
at a neutrino factory with near and far detectors that can identify the lepton charge, and point out 
observables that can discern the sterile mixing signals. We find that clean identification of sterile 
mixing would be possible for ^24^34 > 0.005 and ^14 > 0.06 rad with the current bound of ^13 < 0.2 
rad; a better ^13 bound would allow probing smaller values of sterile mixing. We also generalize 
the formalism for any number of sterile neutrinos, and demonstrate that only certain combinations 
of sterile mixing parameters are relevant irrespective of the number of sterile neutrinos. This also 
leads to a stringent test of the scenario with multiple sterile neutrinos that currently is able to 
describe all the data from the short baseline experiments, including LSND and MiniBOONE. 
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I. INTRODUCTION 



In the framework of the Standard model (SM) of particle physics, there are only three 
neutrinos, one each with the electron, muon and tau flavor. The LEP experiments have 
determined the number of light neutrinos that couple with the Z boson through electroweak 
interactions to be 2.984±0.008 [1], thus closing the door on any more generations of "active" 
neutrinos. However, there still may exist sterile neutrinos that do not have electroweak 
interactions. Though they cannot be detected in the Z decay, they may mix with the active 
neutrinos and hence participate in neutrino oscillations. The extent of sterile neutrino 
participation in solar neutrino data is severely restricted from the neutral current data from 
SNO [2]. The atmospheric neutrino data show that the major contribution to the muon 
neutrino disappearance has to be from i/^ ^ oscillations, however a small admixture of 
sterile neutrinos cannot be ruled out [3, 4]. Short baseline experiments sensitive to sterile 
neutrinos in the ~ 1 eV range [5-8] have given strong upper bounds on sterile mixing. The 
CNGS experiment expects to further restrict the sterile mixing parameter space through 
the i/e,r channels [9]. The MiniBOONE experiment [10] has virtually ruled out any 

effect of sterile neutrinos in the LSND parameter space [11] if there were only one sterile 
neutrino species. However, recently it has been pointed out [12] that with two or more sterile 
neutrinos, it is possible to be in agreement with all the data. 

Even if the LSND results are ignored, so that there is no longer any need for sterile 
neutrinos for explaining the neutrino oscillation data, sterile neutrinos that obey all the 
constraints from the terrestrial experiments can still play a crucial role in astrophysics and 
cosmology [13]. The matter enhanced active-sterile neutrino transformation can have a great 
effect on r-process nucleosynthesis in the core-collapse supernovae [14], and can also influence 
the explosion dynamics [15]. Moreover, the anisotropy inside the exploding supernova may 
be transported outside efficiently by sterile neutrinos, thus helping to explain the large 
observed velocities of pulsars [16]. Sterile neutrinos of mass ~ keV are also excellent dark 
matter candidates: i/MSM (SM with three sterile neutrinos) [17] manages to explain masses 
of active neutrinos, baryon asymmetry of the universe and the abundance of dark matter 
together. The Chandra blank sky observations also allow keV neutrinos to be viable dark 
matter candidates [18]. Such heavy dark matter also helps in the production of supermassive 
black holes [19]. Sterile neutrinos may leave their imprints in the supernova neutrino burst 
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[20, 21], or in the ultrahigh energy neutrino signals observed at the neutrino telescopes [22]. 

The main requirements for the astrophysically and cosmologically relevant sterile neu- 
trinos are thus that they be heavy (m ~ 1-10 eV for r-process nucleosynthesis, and m ~ 
keV for the dark matter candidates) and that they mix weakly with the electron and muon 
neutrino (in order to satisfy the MiniBOONE constraints). In this article, we consider the 
case with one such sterile neutrino, which may influence neutrino oscillation experiments. 
We perform the complete 4-i/ analysis, taking into account all the three additional mixing 
angles of the sterile neutrino Vs with the active ones, and the two additional CP violating 
phases. We only concentrate on heavy neutrinos, such that if 777.4 is the mass of the neutrino 
eigenstate with a dominant sterile component, Am^^ = \m1 — mf\ > 0.1 eV^ for all other 
neutrino mass eigenstatcs z/j. The oscillations due to Amg^ are rather rapid, and can be 
taken to be averaged out in the long baseline data. As a result, the data are expected to be 
insensitive to the exact value of Amg^. However, the additional mixing angles may leave 
their signatures in the data. 

We treat the effects of the sterile neutrino as a perturbation parametrized by a small 
auxiliary parameter A = 0.2. To this end, we represent the active-sterile mixing angles 
^14, 6'24, ^34, the deviation of ^23 from maximality, the reactor angle ^^13 as well as the ratio 
AjTig/Am^^-jn, formally as some power of A times 0{1) numbers, so that a systematic ex- 
pansion in powers of A may be carried out. Averaging out the fast oscillations due to Amg^ 
allows us to obtain simple analytic approximations for the flavor conversion probabilities of 
neutrinos. The expressions thus obtained describe the dependence of relevant conversion or 
survival probabilities on the parameters in a transparent manner. 

We analyze, using analytical as well as numerical means, how the parameters involving 
sterile neutrinos - constrained by the data from solar, atmospheric, and short baseline 
experiments - affect the results at the long baseline experiments. We illustrate this effect 
quantitatively in the case of a neutrino factory setup involving a near and a far detector that 
are capable of lepton charge identification. In particular, we consider the CP asymmetry in 
IX and T channels as the observables and calculate how far the limits on the sterile mixing 
parameters can be brought down. We also consider the electron channel, where signals of 
sterile neutrino mixing can still be established by the counting of the total number of events 
above a threshold. 

The paper is organized as follows. In Sec. II, we explain our formalism of a systematic 
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expansion of all quantities in an auxiliary small parameter A and the use of perturbation 
theory to obtain the neutrino flavor conversion probabilities. In Sec. Ill we examine some of 
the possible signatures of sterile neutrino mixing on the signals at a neutrino factory setup 
with near and far detectors, where we also estimate bounds that can be obtained at such 
long baseline experiments. In Sec. IV, we generalize our formalism to any number of sterile 
species and point out that only certain combinations of the sterile mixing parameters are 
relevant, independent of the number of sterile species. Sec. V concludes. 



II. ANALYTIC COMPUTATION OF NEUTRINO FLAVOR CONVERSION 
PROBABILITIES 



We work in the A-i/ framework, where (z/g, z/^, z/^-, i/g) form the basis of neutrino flavor 
eigenstates and {ui, ^2, t's, i'a) form the basis of neutrino mass eigenstates. The mass eigen- 
states are numbered according to the convention |Am|2| ^ 1^^321 ^ ^^21 > 0> where 
Am?- = — m|. We have Amg^ |Am42|, Am^^m ~ |Am|2| and Am| ^ Aml^. Note that 
the sign of Am32 is as yet unknown, a positive (negative) Amgg corresponds to the normal 
(inverted) mass ordering of neutrinos. 

The mass and flavor eigenstates of neutrinos are connected through a unitary matrix U, 
such that 

i^a = i^ail^i , (1) 

where a e {e, 11, r, s} and i e {1, 2, 3, 4}. The mixing matrix U may be parametrized as 

U = Uuieu, Su) Uu{eu, O) C/24(^24, C/23(^23, 0) C/l3(^13, <^13) C/l2(^12, 0) , (2) 

where Uij{9ij, 5ij) is the complex rotation matrix in the i-j plane, whose elements [C/jjlp^ are 
defined as 

cos 9 p — q — iorp — q — j 

1 P — Q ^ P — Q j 

sin 9e~^^ p = i and q = j (3) 



sin ^e*^ P = j and q 



otherwise . 

The limit when the sterile neutrino is completely decoupled - or when it does not exist - 
is obtained simply by putting the mixing angles ^14,^24 and ^34 to zero. In this limit, U 
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in (2) reduces to the standard Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mixing 
matrix. Since we are interested only in oscillation experiments, we neglect any Majorana 
phases. 

We expect ^14,^24 and 634, the mixing angles involving the sterile neutrino, to be small. 
Indeed, though the 4z/ analysis of atmospheric neutrinos give a rather weak bound of 624 ~ 
|W^4p < 0.19 [4], short baseline disappearance experiments [5] constrain < 0.013, whereas 
the short baseline appearance experiments [6-8, 10] give a bound of ^14^24 ~ \i^ed^ij,4:\ < 
0.02. The atmospheric neutrino data restrict the deviation of $23 from maximality to be 
< 0.15 rad [23], and the CHOOZ data [24] combined with solar, atmospheric and KamLAND 
experiments constrain 613 to be less than 0.2 rad [25]. In order to keep track of the smallness 
of quantities, we introduce an auxiliary number A = 0.2 and define the small parameters to 
be of the form aA". This allows us to perform a systematic expansion in powers of A. For 
the sterile mixing angles, we define 

^14 = X14A ^24 = X24A 9s4 = X34A , (4) 
whereas for the active mixing angles, we define 

TT ~ TT 

^13 = X13A , ^23 = ^ + ^23 = ^ + X23A . (5) 

Here, all the Xij ^^e taken to be 0{1) quantities. We also treat the solar mixing angle, 
9-12 ~ 0.6, as an 0{1) quantity. The limits on the other ^^s mentioned above translate to 
X24 < 0.6, X14X24 < 0.5, X23 < 0.75 and Xis < 1- 

In the long baseline neutrino experiments, the dominating term in flavor conversions 
oscillates as sm'^[Aml^j^L/{4:E)]. Owing to the small value of AmQL/(4£'), the oscillations 
due to AuIq do not have enough time to develop, and the effect of AuIq may be viewed as a 
perturbation to the dominating Aml^^ oscillations. We treat the ratio AmQ/\Aml^^\ ^ 0.03 
as a small parameter, and define 

AmlJAml2 = CA' ■ (6) 

Note that ( is positive (negative) for the normal (inverted) neutrino mass ordering. 

When neutrinos pass through the earth matter, there are matter effects that give rise 
to an effective potential T4 — V^GpNe for the electron neutrino as compared to the other 
neutrinos by virtue of the its charged current forward scattering interactions. Here Gp is 
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the Fermi constant and is the number density of electrons. In addition, all the active 
neutrinos also get an effective potential Vn — —GpNn/V^ compared to the sterile neutrino 
by virtue of their neutral current forward scattering reactions. Here A^^^ is the number 

density of neutrons. For antineutrinos, the signs of Vg and Vn are reversed. The effective 
Hamiltonian in the flavor basis is then 

/-Amii \ 


\ml^ 
V Aml^J V 
where Ae(n) = 2i?K(n), and Uq is the mixing matrix in vacuum, whose form is given in (2). 
Let Hf be diagonalized by a unitary matrix U^, such that 

HD=UlHfUm, (8) 

where Hd is the diagonal matrix. The elements [HdIh, being the eigenvalues of Hf, give the 
relative values of mf/{2E), where rfii are the effective masses of the interaction eigenstates 
in matter. If we assume that the density encountered by the neutrinos during their passage 
through the earth is a constant, the flavor conversion probabilities may be written in terms 
of mf and the elements of Um as 



Pa/3 = P{l^a 



X][Wm]m[Wm]^jexp 



2E 



(9) 



This approximation is valid as long as the neutrino trajectories do not pass through the 
core, and the neutrino energy is not close to the 613 resonance energy in the earth. 

In order to calculate Um, it is convenient to work in the basis of neutrino mass eigenstates 
in vacuum. The effective Hamiltonian in this basis is 



Hv = ul Hf Uq 

/-Amii 



1 

2E 



LV 


















Am^2 



\ 



Amis/ 



( Ag + An Q \ 









A„ 
A„ 
0/ 



(10) 



which can be diagonalized by the unitary matrix U defined through 



Urr,. = UM 



(11) 



such that 

Hd = C/^ Ul HfUoU ^U^ H^U . (12) 

Using the formal representation of the elements of Uq as well as Am^i in terms of A as shown 
in eqs. (4), (5), and (6), the matrix Hy can now be expanded formally in powers of A as 

Hv = ^ [ho + Xh + X'h + 0{X')] . (13) 

The elements of /io,i,2 are functions of all the neutrino mixing angles, mass squared differences 
and CP violating phases in general; the exact expressions are given in Appendix A. All the 
elements of the matrices hi and /12 are of 0{1) or smaller, so that the techniques of time 
independent perturbation theory can be used to calculate the eigenvalues and eigenvectors 
of Hy that are accurate up to 0{X'^). The complete set of four normahzed eigenvectors gives 
the unitary matrix U that diagonahzes Hy through eq. (12). Using eq. (11), one obtains the 
unitary matrix Um that diagonahzes Hf through eq. (8). The matrix Um and the eigenvalues 
of (or Hf) allow us to calculate the neutrino flavor conversion probabilities from eq. (9). 

The flavor conversion probabilities of neutrinos, accurate to (9(A^), obtained by assum- 
ing the neutrinos to travel through a constant matter density, are given in Appendix A. 
These expressions seem rather complicated. However, we can make certain approxima- 
tions that will simplify these expressions and bring forth some important physical insights. 
Since we are interested in heavy sterile neutrinos, we may take |Am|2| <^ |Am|2|- Also, 
since \Aml2L/E\ ~ 0{1), we have \Aml2L/E\ S> 1 and the oscillating terms of the form 
cos{Aml2L/E) may be averaged out. In the long baseline experiments, we are interested 
in the energy range 1-50 GcV. Even at the higher end of the energy spectrum, taking the 
density of the earth mantle to be w 5 g/cc, we get Ae 2 x 10~^ eV^ and A„ fti — 1 x 10~^ 
eV^ for neutrinos, so we also approximate |^e,n| <^ I ^"^42 1 wherever appropriate. With 
these approximations, the neutrino flavor conversion (or survival) probabilities for an initial 
ly^ may be written as 

P^^2gj,Al /"l^^:l^^ +0(.X^), (14) 

Pnn ~ cos^ A32 + 4^^3 sin^ A32 - A21 sin^ 6*12 sin 2 A32 

+7A-^Sr-N2 {-2A32 COS A32 Sin Ae sin(Ae - A32) + A^iA^ - A32) sin 2A32} 

(Ae - A32)^ 

-29l^ cos^ A32 + 2e24^34A„ COS ^24 siu 2A32 + C>(A^) , (15) 
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p, 



sm' 




'32 



(^24 + ^34) sin^ ^32 - ^24^34 (2A„ COS 824 + sin 824) sin 2A32 + 0{\^) , 



(16) 



where we have defined the dimensionless quantities A^j = Am^jL/^AE) and Ae,„ = 
Ae,nL/{4:E) for convenience. The following observations may be made from the above ex- 
pressions: 

• The leading 0{1) terms are of the form sin^ A32 or cos^A32, corresponding to the 
dominating atmospheric neutrino oscillations. There is no subleading term of (9(A). 

• For P^e, there is no sterile contribution up to C(A^). Indeed, the leading order sterile 
contribution to P^e is proportional to ^14^34, which is 0{X'^). 

• In the expression for P^^ or P^^, the first line contains the leading oscillating term 
as well as the subleading terms due to the deviation of 6*23 from maximality and due 
to the nonzero value of Am2i- The next line gives the contribution from 9^^, which 
matches the one obtained in [26] . The last line contains the contribution from sterile 
neutrinos. Whereas it is enough to have either ^24 or ^34 nonzero for the sterile mixing 
to have an effect on P^^, the sterile contribution to P^^ will be present only for nonzero 



• Only one CP violating phase, ^24, is relevant for the flavor conversion probabilities up 

to this order. The phases ^13 and ^14 appear only at 0{X^) or higher. In particular, 
the CP violating terms proportional to (Am2i/Am|2)6'i3, as given in [26], are absent 
since they are of O(A^). 

• Note that the leading sterile contribution at the long baseline experiments is found 
to be at (9(A^). This may be compared with the CP violation in the active sector, 
whose leading contribution appears at 0{X^) and the short baseline appearance ex- 
periments, whose positive results would appear only at (9(A^) or higher. The O(A^) 
sterile contribution to P^t, which is proportional to sinA32, is absent in the short 
baseline appearance experiments where in general IA42I ~ and IA32I <^ 1, so 
that sin A32 ~ 0. 
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• When Ae ~ A32, the 613 contribution is enhanced due to the factor (Ag — A32)~^. The 
analytical approximation is expected to fail in this region since even the higher order 
terms in ^13 may become significant. 

• The analytic expressions are not expected to be valid for large L/E where A21 would 
become 0{X) and higher order terms in A21 would also contribute to the probability 
in (15) and (16) at ^(A^). 

• The probabilities in (14), (15) and (16) do not involve Am^^, and have no information 
on whether the mainly sterile neutrino 1^4 is heavier or lighter than the other three. 
This is due to our approximation of averaging out the fast oscillations due to Am^^. 
This approximation will be more and more accurate as Am^^ increases. 



e.n 



The probabilities for the antiparticles are obtained simply by replacing „ —A 
and 5ij — > — The sterile contribution to the CP violation is therefore given by 

P^,^, - ~ (i^MM - Pmhv + 4^24^34A„ cos 824 sin 2A32 , (17) 
Put - Pfif ~ {Put - Pjifjiv - 46*24 6*34 A„ cos ^24 sin 2A32 - 26'246'34 sin ^24 sin 2A32 . (18) 

The CP violating contribution of sterile neutrinos to P^^ — Pfifi is entirely from the earth 
matter effects, whereas for P^^- — P^^, the contribution comes from both the earth matter 
effects (through the A„ term) as well as the vacuum mixing matrix Uo (from the sin524 
term) . 

For an initial Vg, the relevant neutrino flavor conversion probabilities are 

1^ 

(Ae - A32)2 

P.,^2el,Al,'^^^^^^ + OiX^), (20) 

Pe. « 2el,Al /'f^^^-^^^^ + 0(A3) , (21) 

where we have used the approximations |Ae,n| -C IA42I, and have averaged out terms that 
oscillate as fast as sinA42. The complete expressions accurate to 0{X^) may be found in 
Appendix A. Clearly, sterile neutrinos have no effect at this order on these probabilities 
except on Pee, and there is no sterile contribution to the CP violation in any of these three 
channels. 



Pee^l- AOfAl^ '^y- . ^.f ^ - 2^?4 + 0{X') , (19) 
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FIG. 1: Probabilities Pun and Pur as functions of energy, and the comparisons with the analytic 



expressions in eqs. (15) and (16). In all the plots, we take Attiq = 8 x 10 ^ eV^, A 



atm 



2.5x10-3 



eV2, 6*23 = 45°, and 6*12 = 33.2°. The magenta (dotted) curve corresponds to the situation with 
no sterile contribution and vanishing ^13. The blue (red) curve with rapid (extremely rapid) 
oscillations corresponds to Amgi- = 0.1 (1.0) eV^, with On = 624, = O34 = 0.2 rad and ^13 = 0. The 
black curve that passes through the rapidly oscillating curves denotes the analytical approximation, 
which is independent of the value of Arrig^ since the high frequency oscillations are averaged out. 



The green (dashed) curve represents the situation with 614 = O24 = ^34 = 0, but ^13 = 0.2 rad. 
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We demonstrate the validity (and limitations) of our analytic approximations in Fig. 1, 
where we show P^n and P^r as a function of energy for three baselines, 3000 km, 7000 km 
and 10000 km. In each panel, we show the probabilities with $14 — 62A — ^34 = 0.2 rad 
and ^^13 = 0, for Am^^. = 0.1 cV^ and Am^j = 1 eV^: the complete 4- neutrino numerical 
simulation with the Preliminary Reference Earth Model (PREM) [27] for the density of the 
earth, as well as our analytical approximation that uses the average density along the path 
of the neutrino and averages out the high frequency approximations.^ In order to estimate 
whether nonzero ^13 can mimic the signatures of sterile mixing, we also show the probability 
for all the sterile mixing angles vanishing, but ^13 = 0.2 rad. 

The following observations may be made: 

• The analytical approximation agrees well with the average of the exact numerical 
results for L = 3000 km and 7000 km. For L = 10000 km, though the analytic 

approximation predicts the qualitative behavior of the averaged probabilities, the exact 
numerical values have an error of ~ 5%. This is due to the large L making A21 ~ 
so that higher order terms in A21 contribute to the probabilities (15) and (16). 

• The dominant effect of the sterile contribution is to pull down the value of P^^, which 
mimics the deviation of ^23 from its maximal value. Such a mimicking is also possible 
through a nonzero ^13, however the effect of may be significantly larger, beyond 
what is possible with the current limit on ^13. Moreover, at energies much larger than 
the ^13 resonance, the ^13 contribution is suppressed by the factor A32/(Ae — A32) 
in earth matter, whereas the sterile contribution does not undergo any suppression 
since |A„| <^ IA42I in the whole energy range of interest. One therefore expects that 
distinguishing the sterile contribution would be easier at high energies. 

• Sterile contribution to P^^ as well as P^r is larger at longer baselines, due to the A„ 

term present in (15) and (16), which increases with increasing L. On the other hand, 
at low L/E values, the sterile contribution to P^,- is highly suppressed by the factor 
sinA32 in (16). 

^ For a baseline of 10000 km, we only show Amg^ = 0.1 eV^, otherwise the oscillation frequency would be 
too high. 
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III. SIGNATURES AT LONG BASELINE EXPERIMENTS 



The analytical expressions (14)-(16) indicate that at > 10 GeV where |Ae| » IA32I, 
the contribution of the currently unknown ^13 is suppressed by a factor ~ A32/Ae. There 
is no such suppression for the sterile contribution, since |Ae_„| <C IA42I for £■ < 50 GeV. 
For E ~ 5-10 GeV, the earth matter effects cause an enhancement of ^13 through the 
factor A32/(Ae — A32). This energy range is therefore unsuitable for searching for a sterile 
contribution to the conversion probabilities. At £^ < 5 GeV also, since the contribution due 
to the currently unknown ^13 is at least of the same order as the maximum allowed sterile 
contribution, discriminating between ^13 and sterile contributions to the probabilities would 
need data from more than one experiment. A high energy neutrino experiment is therefore 
preferred. 

In order to demonstrate the capability of future long baseline experiments in distinguish- 
ing the sterile neutrino contribution to the neutrino flavor conversion probabilities, we choose 
a typical neutrino factory setup [28], with a 50 GeV muon beam directed to a 0.5 kt "near" 
detector 1 km away, and a 50 kt "far" detector 7000 km away. The detectors may be mag- 
netized iron calorimeters [29] , which can identify the charge of the lepton produced from the 
charged current interaction of the neutrino or antineutrino. The number of useful muons in 
the storage ring is taken to be 1.066 ■ 10^^, which corresponds to approximately two years 
of running with /i~ and /x"*" each at the neutrino factory, using the NuFact-Il parameters in 
[30] . We implement the propagation of the neutrinos through the earth using the 5-density 
model of the Earth, where the density of each layer has been taken to be the average of 
the densities encountered by the neutrinos along their path in that layer with the PREM 
profile [27]. We take care of the detector characteristics using the General Long Baseline 
Experiment Simulator (GLoBES) [31]. This includes an energy resolution oi (Te/E = 15%, 
an overall detection efficiency of 75% for all charged leptons, as well as additional energy 
dependent post-efficiencies that are taken care of bin-by-bin. We assume perfect lepton 
charge identification, and neglect any error due to wrong sign leptons produced from the 
oscillations of the antiparticles. These can be taken care of in the complete simulation of 
the detector once its detailed characteristics are known. 
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FIG. 2: The asymmetries A^{E) and Ar^E) as functions of energy at a neutrino factory. The 
band corresponds to ahowed values of the asymmetries without any sterile mixing, with ^23 1 ^13 
and 613 allowed to vary over all their allowed ranges, and with both the normal (NH) and inverted 
(IH) hierarchies. The plots for showing the dependence on sterile components are with 623 = vr/4, 
^24 = 0, and Am42 = 0.1 eV^. The results will not change if Am\2 has higher values. No significant 
dependence on ^14 is expected from (15) and (16), hence we use On = 0. The errors shown are 
only statistical. 



In Fig. 2, we display the asymmetries 

where Ni (Ni) is the number of (£+) observed at the near or far detector. These asymme- 
tries roughly correspond to Af_i ~ P^^ — P^^ and Ar ~ Prr — Pff, where the events observed 
in the near detector act as a normalizing factor, and help in canceling out the systematic 
errors due to fluxes, cross sections and efficiencies in each energy bin. Note that we do not 
expect any at the near detector, hence the number of events of at the far detector 
needs to be normalized to the number of events of fi^ at the near detector. 

In the absence of any sterile neutrinos, and in the limit of vanishing 6*13, the asymmetries 
A^ and At vanish, as can be seen from (15) and (16). The ^13 contribution is indeed 
suppressed at high energies, as discussed above. In the figure, we show a band corresponding 
to the possible signals in the absence of any sterile neutrinos, where we vary over the allowed 
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values of the angles 623, 613, the CP phase ^13 and both the normal as well as inverted 
mass ordering. For Aml^^, AtUq and Ou we only take the current best-fit values, since the 
variation in these parameters is not expected to cause any significant change in our results. 
We choose to take 624, = O34 and ^24 = for illustration, since from (17) and (18) we expect 
the asymmetries to be identical in magnitude and proportional to the product ^24^34 with 
vanishing ^24. Any discrepancy between these two asymmetries would indicate a nonzero 
^24, and hence CP violation in the sterile sector. The third sterile mixing angle, 9u, is taken 
to be vanishing since it is not expected to affect the relevant neutrino conversions. 

It may be observed from Fig. 2 that for E > 15 GeV, the sterile contribution results in 
an deficit (excess) of the asymmetry for normal (inverted) hierarchy in the /i channel. In 
the r channel, the situation is the reverse. This is as expected from our analytic expressions 
(17) and (18). The asymmetry integrated over energy may therefore be expected to serve 
as an efficient discriminator between the scenarios with and without sterile neutrinos. In 
Fig. 3, we show the integrated asymmetries 

~ _ Nl^'{E > 15GeV) N^^\E > 15GeV) 



iV--(E > 15GeV) N^'^'^iE > 15GeV) ' 

~ ^ Nl^^jE > 15GeV) _ nI^\E > 15GcV) 

7V°«ar(£; > i5GeV) ]V;'^'(£; > 15GeV) ' ^ ^ 

The figure indicates that for ^24^34 > 0.005, the sterile contribution to neutrino conversions 
can be discernable from the three neutrino mixing results. The width of the band is de- 
termined essentially by the allowed range of ^13. If the value of ^13 is bounded further, 
the reach of neutrino factories for the sterile mixing is enhanced. In addition, the actual 
value of ^13 also affects the discovery potential of sterile mixing by infiuencing the integrated 
asymmetries A/^, At, as shown in the figure. Note that since the asymmetries depend on 
the sign of Amgg, sterile mixing also makes it possible to distinguish between normal and 
inverted hierarchies. 

If we have a 50 kt detector that can detect e~ je^ and identify their charge^, we can use 



^ Charge identification is needed in order to get rid of the error due to misidentification of the wrong sign 
leptons produced due to v^^ — > or i/^ — > oscillations. A magnetized iron calorimeter with thin iron 
strips, or a liquid Ar detector [32], may serve the purpose. If charge identification is not possible, as in a 
water Chcrenkov detector for example, the background due to the wrong sign lepton will have to be taken 
into account. 
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FIG. 3: The integrated asymmetries and At as functions of sterile mixing parameters at a 
neutrino factory. We use 624 = 034 and 824 = 0. The rest of the parameters are the same as in 
Fig. 2. The statistical errors are smaller than the circles shown in the plots. 



the observable 



and the integrated quantity 



yynear 
e 



Nl^'iE > 25GeV) 



7^, 



(24) 



(25) 



JYnear(^ > 25GeV) 

for detecting the sterile neutrino contribution. Note that there is no difference between the 
two hierarchies, or between z/g and z/g, as far as the expected probabilities are concerned. 
From Fig. 4, it may be seen that for 614 > 0.06, the sterile mixing signals can be clearly 
discerned. If the bound on ^13 becomes stronger, even smaller values of 614 may be identified. 
On the other hand, an higher actual value of 6*13 helps in the identification of sterile mixing 
even at lower 6*14 values. 

The "platinum" channel P^e at the neutrino factories is not affected by the sterile mixing, 
not just to O(A^), but even at 0{X^). Indeed, going to one higher order in the A-perturbation, 
we get 

A^^ sin2(Ae- A3 



P, 



^32 J 



fie 



26*^3 



(Ag - A32)2 

-2^i3A2iA32 sin (2^12) sin A32 



4^13(^^23 



-^13 

COS (Ae - A32 
(Ae-A 



'13 j 



sin Af 



32 j 



(26) 
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FIG. 4: The observables TZe{E) and TZe at a neutrino factory, where e~ /e~^ and their charge may 
be identified. The active neutrino mixing parameters are the same as that in Fig. 2. The sterile 
mixing parameters are taken to be 624 = ^34 = and Am42 = 0.1 eV^. Any increase in Attt-Ij, 
or nonzero value of ^24/^34 are not expected to have any significant effect on this observable. The 
result is insensitive to sgn{Am'^2)- 



For getting P^g, one just needs to replace Ag — Ag and 6 is —Si^. This channel is 
therefore not expected to be useful in putting constraints on sterile mixing. On the other 
hand, it is free of any sterile contamination to 0{X^), and is therefore suitable for determining 
the parameters in the standard three flavor analysis. 



IV. GENERALIZATION TO ANY NUMBER OF STERILE NEUTRINOS 

If the LSND results [11] are taken to be valid, a single sterile neutrino is not enough to 
describe all the data from short baseline experiments. However, two or more sterile neutrinos 
with Amj^ ~ 1 eV (j > 3) and |WejW^j| ~ 0(0.01-0.1) are consistent with all data [12]. 
Some avenues for probing the mixing parameters and distinguishing between different mass 
orderings in such a case have already been suggested [21, 33]. It is therefore desirable to 
extend our formalism to more sterile neutrinos. 

The analytical treatment in Sec. II for the case of one sterile neutrino may be generalized 
easily to any arbitrary number n of sterile neutrinos. The (3 + n) x (3 + n) mixing matrix 
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U may be written in the block form as 



U 



nxn 



[W 



SAl 



[03X3 [0j3x„ 



nx3 



[W 



SS] 



(27) 

where = C/23(^23,0) C/i3(^i3, (^13) ^^12(^12,0) is the standard mixing matrix for three 
active neutrino flavors, and V^^ is the matrix that mixes the n sterile neutrinos among 
themselves. Since the assignment of "flavor" eigenstates to the sterile species is arbitrary, 
we choose the basis such that the flavor and mass eigenstates of the sterile neutrinos coincide 
in the absence of any activc-stcrilc mixing, i.e. V^^ = Inxn- The matrix W parametrizes the 
mixing between active and sterile states, and in general may be represented by a product of 
matrices Uij{9ij, Sij) as defined in eq. (3), with i < 3 and j > 3. 

In addition, we assume that all the active-sterile mixing is small, which is borne out by the 
recent 3+2 neutrino fit to LSND, MiniBOONE as well as the short basehne disappearance 
data [12]. This allows us to write 



( 



\ 



I 



= A 



[X'''^]lxn 



(28) 



If terms of C(A^) and smaller are neglected, the unitary matrix W may be written in its 
most general form as 



J3x3 



x3 



3xn 



+ C(A2 



(29) 



The net leptonic mixing matrix U in (27) can then be written as 



U 



2 



V 



AA 



3x3 



nx3 



(30) 



For the active mixing angles in V^^^ we use the same A-expansion as in eq. (5), i.e. ^13 = X13A 
and ^23 = 7r/4 + ^23 = 7r/4 + X23A. We also treat Aml^/Amga to be a small quantity, and 
denote it formally by Am|i/Am|2 = C,y?-, as in eq. (6). The quantities C^Xij ^ well as all 
the elements of X^^ are taken to be Oil) parameters. 
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Following the same systematic expansion procedure delineated in Sec. II in the case of 
one sterile neutrino, we obtain the neutrino flavor conversion (or survival) probabilities for 
an initial i/^ beam to be 

P,«20iSg^^ + 0(A=), (31) 



P^^ fti cos^ A32 + 46*23 sin^ A32 - A21 sin^ 6*12 sin 2 A32 

+ f ,2 {-2A32 COS A32 sin Ae sin( Ae - A32) + Ae( A^ - A32) sin 2A32} 

(Ae - A32)^ 

-2[W^^{W^'^)^] cos^ A32 + 2Re[M/"^(iy^^)^] A„ sin 2A32 + C(A^) , (32) 



P^,r ^ sin^ A32 - 46'23 sin^ A32 + A21 sin^ 612 sin 2 A32 

^3 {2A32SinA32COsAeSin(Ae- A32) - Ae(Ae- A32)sin2A32} 



(Ae - A32)2 

- {[W''^{W^'^)^ + [W^^{W^^)^]) sin^ A32 

-2Re[W^^^(W^''^)^] A„ sin 2A32 - Im[W^^^(W^'^^)^] sin 2A32 + ^(A^) . (33) 

Here we have assumed |Am|2|, |^e,n| <^ lAm^gl, and have averaged out the oscillating terms 
of the form cos(Am|2-^/-E), as before. The sterile contribution to the CP violation in these 
channels is then 

- Pp-^ ~ (Pm^ - Paa)3. + 4Re[W^-^(W^'^^)t] A„ sin 2A32 , (34) 
P^r - Put ~ {Pi^r - Pfifh. - 4Re[iy^^(iy'^^)t] A„ sin 2A32 

- 2lm[W^^{W^^)^] sin 2A32 . (35) 

For an initial beam, the corresponding flavor conversion probabilities are 

Pee 1 - AelA ^'f^^'-^'^ - m^'iW^')^] + 0{X') , (36) 
[lAe — A32j 

Pe, - 2^i3Al2^^^^^^^ + 0{X^) , (37) 

Pe. ^ 2el, ^'f}^\^^ + 0{X') . (38) 



The mixing matrix 14 in (30) reduces to the 4x4 mixing matrix U (2) in the case of one 
sterile neutrino simply by taking n = 1 and using the substitution 

W^^ eue~'^'^ , W^"^ 0246-'^^^ , W^^ Ou ■ (39) 
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As a result, the bounds obtained on 614, 624, O34 and 824 in the 4-neutrino analysis can be 
directly translated to bounds on the combinations [W^iW^Y], [W^'^(Wi'^)% [W^^(W^^y] 
as well as the real and imaginary parts of Note that the expressions (14)-(21) 

obtained in the special case of only one sterile neutrino can be obtained from the general 
expressions (31)-(38) simply with the substitutions (39). Specifically, the bounds obtained 
on ^24^34 in Sec. Ill using the observables Afj,,AT are simply bounds on Re[W^'''^(W^^'^)^]. 
Similarly, the bound obtained on Ou through is simply the bound on [VF'^'^(iy^'^)^]^/^. 

The above argument also implies that, at least in the region of validity of our ana- 
lytic approximations, the only combinations of active-sterile mixing parameters that may 
be bounded by data arc the four quantities [W-'^iW^'^y], [W^iW^Y], [W^'^{W^^y] and 
[W^^'^(iy^'^)^], irrespective of the number of sterile species. For example, in the 3-1-2 scenario, 
the mixing matrix U may written as 

U — U45{945, ^45) • U35{935, ^35) • 1/25(025, S25) ■ Ui5{9i5, 615) ■ t/34(6'34, ^34) " 

C/24(^24, 524) ■ UuiOu, Su) ■ C/23(^23, <^23) ' C/l3(^13, <^13) ' C/l2(^12, <^12) , (40) 

where 645 — 0, and 9ij ~ 0{X) for j > 3. One may, in addition, choose some of the phases 
5ij to be vanishing by proper redefinitions of leptonic phases. With the mixing matrix U in 
(40), the substitution 

-iS25 



\ 



6246-'^'' 0256 

'34e-'^'' 635 



(41) 



would give the relevant combinations of the sterile mixing parameters: 



[W^^iW''^)^ 



15 ' 



On A + 9^ 



25 ' 



^34 + ^35 > 



9249346'^^''-^''^ + 925935e'^^''-^''^ 



(42) 



The expected bounds obtained in Sec. Ill then would correspond to 

^24^34 C0S(524 - ^34) + ^25^35 cos(525 " S35) < 0.005 , ^J 9^ + 9^, < 0.06 . (43) 

These bounds will act as a stringent test of the scenario with multiple sterile neutrinos [12] . 
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V. CONCLUSIONS 



Heavy sterile neutrinos may play an important role in astrophysics and cosmology, for 
example in r-process nucleosynthesis or as dark matter. Neutrino oscillation experiments, 
mainly the short baseline ones, have already put severe constraints on the extent of mixing 
of these sterile neutrinos with the active ones. If the LSND results are taken to be valid, at 
least two sterile neutrinos are in fact needed to describe all data. 

Our aim in this paper is to check whether the sterile neutrinos so constrained can still 
give rise to observable signals at future experiments, and whether these signals can be 
cleanly identified in spite of our current lack of knowledge of all parameters in the mixing 
of three active neutrinos. This would lead to an estimation of bounds on the sterile mixing 
parameters that can be obtained with neutrino oscillation experiments. 

The number of neutrino mixing parameters increase quadratically with the number of 
neutrinos, and only certain combinations are expected to be relevant for neutrino flavor 
conversions. In order to identify these combinations in an analytically tractable manner, 
we exploit the smallness of certain parameters to carry out a systematic expansion in an 
arbitrarily defined small parameter, A = 0.2. The small quantities ^14, ^24, ^34, ^13, ^23 — 7r/4, 
and Am0/Am^tjjj are formally written as powers of A times 0{1) numbers, and neutrino 
conversion probabilities correct to (9(A^) are obtained using techniques of time independent 
perturbation theory. We also neglect terms proportional to Am1^^^/Aml2, and average away 
the fast oscillating terms like coslAm^L / E) since \Aml2L/E\ 1 in typical long baseline 
experiments. 

It is observed that the conversion probabilities P^g) Pen or Per get no sterile contribution to 
C(A^). For P^^ and P^r, sterile mixing gives contributions proportional to and {Ol^ + O^^) 
respectively. In addition, there is a CP violating contribution proportional to ^24^34 to both 
these quantities. The survival probability P^e gets modified simply by a term proportional to 
9I4. There is no dependence on the mass of the sterile neutrino, since all the terms containing 
Am|2 are averaged out. It is observed that as long as the neutrinos do not pass through 
the core of the earth, the probabilities obtained through our analytic approximations match 
the exact numerical ones rather well. Note that the sterile contribution to the conversion 
probabilities at long baseline experiments appears at (9(A^), which is at a lower order than 
the appearance of CP violation in the active sector or the sterile contribution to short 
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baseline appearance experiments. 

Whereas the contribution due to the currently unknown ^13 decreases at high energies due 
to the earth matter effects, the sterile contribution stays almost constant, and therefore the 
energy range E = 10-50 GeV is suitable for distinguishing the sterile "signal" above the 613 
"background" . The CP violating part of the sterile contribution builds up with increasing L, 
and hence longer baselines are preferable. This naturally leads to the consideration of neu- 
trino factories with = 50 GeV and basehne of a few thousand km as the desirable setup, 
with lepton charge identification capability and a near detector for calibration purposes. 

For illustration we take the far detector to be near the magic baseline of ~ 7000 km, 
and choose three obscrvablcs, and Ar that correspond to the CP asymmetries in the /i 
and r channels respectively, and TZ^, which corresponds to the disappearance in the electron 
channel. The background in these channels is obtained by varying over the unknown values 
of ^13 , 023 and the CP phase S13. It is observed that the signal rises above this background 
for A^ and Ar when ^24^34 > 0.005, and for Tie when Ou > 0.06 rad. The range of 9i4 probed 
is limited mainly by the unknown value of ^13. The limit on ^13 may be brought down by a 
factor of two or more at the reactor experiments like Double CHOOZ [34] or Daya Bay [35], 
and indeed at the neutrino factories themselves [36]. The values of 9i4 that can be probed 
then decrease by approximately the same factor. 

Note that we have only chosen to analyze a few specific observables whose dependence 
on the sterile mixing is analytically transparent. A complete analysis that fits for all the 
parameters simultaneously may give rise to more stringent constraints. The long baseline 
experiments thus have the capability of tightening the limits on the sterile mixing angles 
by almost an order of magnitude over the current ones, or identify sterile neutrinos if their 
mixing is indeed above such a value. Note that if the sterile mixing is identified through Af^ 
or At, the neutrino mass hierarchy - normal vs. inverted - is also identified. 

In the light of the recent results that show that LSND, MiniBOONE and the earlier null- 
result short baseline experiments can be consistent if the number of sterile neutrinos is two or 
more, we have also extended our formalism to include any number of sterile neutrinos. The 
number of distinct combinations of sterile mixing parameters remains the same, irrespective 
of the number of sterile neutrinos. We give explicit expressions for such combinations, and 
the neutrino conversion probabilities in terms of them. The limits obtained on through the 
4-i/ analysis can easily be translated to the corresponding combinations of these parameters 
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in the general case. Indeed, the bounds on the sterile mixing parameters obtained from 
the measurements described in this paper would act as stringent tests of the scenarios with 
multiple sterile neutrinos. 
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APPENDIX A: FLAVOR CONVERSION PROBABILITIES USING PERTUR- 
BATION THEORY TO SECOND ORDER 



In order to calculate the neutrino conversion (survival) probabilities in the presence of a 
sterile neutrino, we define an auxiliary small parameter A = 0.2, write all the small quantities 
as aA" where a and n are some constants, and then perform a formal expansion of the effective 
Hamiltonian in powers of A. This enables us to use the second order perturbation theory to 

get results accurate to (9(A^). 

We have defined the small quantities in the problem as 

^14 = X14A , O24 — X24A , Ou — X34A , (Al) 

^i3 = Xi3A, ^23 - 7r/4 = X23A , Am|/Am2,^ = CA^ (A2) 

As argued in Sec. II, we need to diagonalize the effective Hamiltonian if^, given in eq. (10). 
This Hamiltonian matrix may be expanded in powers of A as 



H„ — 



32 



2E 



[ho + Xhi + X^h2 + 0{X^)] 



(A3) 



Here, the leading term is 

/ a„ + ae cos^ 6*12 ae cos 9i2 sin 612 



ho = 







cos 9i2 sin O12 + sin^ 612 



\ 














a„ + 1 
a / 



(A4) 
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where ae,n = ■^e,n/^nT'l2 = ^"^42/'^"^32 ~ ='=^"^st/'^"^atm- We take the neutrinos to 

be traversing through a constant matter density, so that ae,n are constants. 
The subleading term in (A3) is 



hi 



f aexncoseue-'^'' K + an)xi4e-^'^" cos^ia - ^sin^i2(x24e-^*'* - X34) \ 

aeXis sin Oue-'^'^ (ae + a„)xi4e-''^" sin 612 + ^ cos ei2(x24e-^'^2* - X34) 
. . ^(X24e-*^^^ + X34) 

... 



(A5) 



The matrix hi is hermitian, so wc do not write its lower triangular elements for the sake of 
brevity. Note that all the elements of hi are 0{1). 

The expression for the matrix h2 in (A3) is rather complicated, we just give its ten 
independent elements separately here for the sake of completeness. The diagonal elements 
are 

hl^ = -A32C - [OeXlS + («e + On) XmJcOS^ 612 - y (x24 + X34 " 2X24X34 COS ^24) siu^ 612 

-\/2 (ae + an) Xi4 [X34 COS 6u - X24 cos(-(5i4 + ^24)] sin 6^12 cos 6*12 , 



hf = -[aeXl3 + («e + On) XmIsIH^ 6*12 - y(X24 + Xu " 2X24X34 COS S24) COS^ 012 

+V2 (oe + a„) Xi4[X34 COS Su - X24 cos(-5i4 + S24)] siu COS Ou , 



^2 = -y (X24 + X34 + 2X24X34 COS ^24 j + aeXl3 , 



h^ = On (X24 + X34) + («e + ^n) Xl4^ , 



(A6) 



while the off-diagonal elements are 



hl^ = (-[aeXl3 + («e + an) X14] + y (X24 + X34 - 2X24X34 COS S2S) Sm 012 COS 
(Oe + an) 



+ 



+1 



V2 
. {ae + an 



Xl4[X34 COS(5i4 - X24 COs(-(5i4 + ^24)] COS 26^12 

-Xi4[-X34 sin 5i4 - X24 sin(-5i4 + 524)] , 



^13 ^ - i^^^^''" + ^34) ^ XMe-'-^^^ COS ^12 - an {X24e-'''' - X34) sin 612] , 
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hf = - ^^'"""'''^ ^ ^''^ [V2 (ae + an) Xue-''^' sin 6,2 + {x24e~'''' - X34) 008^12], 
/i^^ = ^(;^24e-'5.4 + ^3^) (-Xi3e-''"cos^i2 + X23sin^i2) , 

^24 ^ --^(x24e-^^^^+X34)(xi3e-''" sin ei2 + X23 003^12) , 

/if = -^X23 (X24e-^^- - X34) + (ae + an) Xi3Xi4e^('--'") . (A7) 

Note that all the elements of /i2 are 0{1) or smaller. The dependence on AttIq appears only 
at this order, and only in the element h^^. 

Using the above formal expansion of the effective Hamiltonian, one can compute the 
eigenvalues and eigenvectors of correct up to C(A^) by using the techniques of time 
independent perturbation theory. The complete set of four normalized eigenvectors gives 
the unitary matrix U that diagonalizes through eq. (12). Using eq. (11), we can then 
compute the unitary matrix Um that diagonalizes Hf through eq. (8). The matrix Um and 
the eigenvalues of Hy (or Hf) allow us to calculate the neutrino flavor conversion probabilities 
from eq. (9). The complete expressions, accurate to O(A^), are given below. 

= 2A^xLAL^^g^^ + OiX^) , (A8) 
P^^ = cos^ A32 + 4X\l^ sin^ A32 - X\ sin^ A32 sin 2 A32 

^2y2 ^ 

+ r A a\2 {-2A32 cos A32 sin Ae sin( Ae - A32) + Ae(Ae - A32) sin 2A32} 

(-Ae + A32)^ 

+A'X24 Ql + A'X34 Q2 + A'X24X34 COS (^24 Q3 + 0{X^) , (A9) 

P^r = sin^ A32 - 4A^X23 sin^ A32 + A^C sin^ ^12 A32 sin 2 A32 

+ ^ A a\2 {-2A32 COS Ae Sin A32 sin ( A32 - Ae) + Ae(-Ae + A32) sin 2A32} 

(-Ae + A32)^ 

+A'(X24 + X34) <34 + A'X24X34(C0S ^24 Q5 + Sm ^24 Qe) + 0{X^) , (AlO) 

where we have defined 

^ 4(A„ + A32-A42)2(-A„ + A42)2 

I - (A„(A32 - 2A42) + 2A42(-A32 + A42))^C0S2A32 

+ (A„A32 - 2(A„ + A32)A42 + 2A42')' COS (2A„ - 2A42) 

+2A2A32(A„- A42)(A„ + A32 — A42) sin 2A32 

-2 (A„A32 - 2(A„ + A32) A42 + 2A42')' sin2(A„ + A32 - A42) } , (All) 
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Q2 = 



32 



2(A„ + A32-A42)2(-A„ + A42)2 

|(A„ - A42)(A„ + A32 - A42) sin2A32 

-2A32 cos A32 sin (A„ - A42) sin(A„ + A32 - A42) | , 



(A12) 



Q3 



An (A„(A32 - 2A42) + 2A42(-A32 + A42)) COS A; 



32 



(A„ + A32-A42)2(-A„ + A42)2 

2(A„ - A42)(A„ + A32 - A42) sin A32 
+A32 [- cos A32 + cos(2A„ + A32 - 2A42)] } , 



X 



{ 



(A13) 



Q4 



8(A„ + A32- A42)2(-A„ + A, 



42; 



■|4A„(A32 - 2A42)(A32 - A42)A42 



-4(A32 - A42)'A^2 - 2A^(A^2 - 2A32A42 + 2Ai^) 



+2 



- 2A„(A32 - 2A42)(A32 - A42)A42 + 2(A32 - A42)'Ai2 



+A^(A^2 - 2A32A42 + 2A^2) 



COS 2 A32 



+ A„A32 sin A;: 



^2 



- 8A„(A„ - A42)(A„ + A32 - A42) COS A. 



32 



-4 [A„(A32 - 2A42) + 2A42(-A32 + A42)] Sin (A32 - 2A42 + 2A„)] } , (A14) 



sin A 



32 



X 



2(A„ + A32-A42)2(-A„ + A42)2 
|A„[A„(A32 - 2A42) + A42(-A32 + A42)] X 

[4(A, - A42)(A„ + A32 - A42) COS A32 + 2A32 sin A32] 



+2 



2 a2 
■42 



- 2A„(A32 - 2A42)(A32 - A42)A42 + 2(A32 - A42)'A 

+Al{Al^ - 2A32A42 + 2A%)^ sin (2A„ + A32 - 2A42)} , 



Qe = 



-4(A32 - A42) A42 sin (A„ - A42) sin A32 sin (A„ + A32 - A42) 



(A„ + A32 - A42)(-A„ + A42) 
Here we have used the shorthand 

AnL 



Ae = — ^ , 



A„ = 



^32 = — 77:; 



_ Am%L 



AE, 



(A15) 



(A16) 



(A17) 
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The probabilities Pea are 

P -1 ^/?2 A 2 sin^(Ae- A32) 4^2 A 2 ^^^^l(^e + A^^^2^ ^ nQ3^ ( A^R\ 
(Ae - Aas)^ (Ae + A„ - A42)^ 



- 2^-^32^^^^^^^ + C!?(A^) , (A19) 
in' (Ae - A3 

(Ae - A32)^ 



Per = 2^^3AL ""l^^\^f^ + 0(A3) . (A20) 



[1] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[2] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 89, 011301 (2002) [arXiv:nucl- 
ex/0204008]; B. C. Chauhan and J. Pulido, JHEP 0412, 040 (2004) [arXiv:hep-ph/0406227]. 

[3] S. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 85, 3999 (2000) 
[arXiv:hep-ex/0009001]. 

[4] G. L. Fogh, E. Lisi and A. Marrone, Phys. Rev. D 64, 093005 (2001) [arXiv:hep-ph/0105139]; 

T. Nakaya [SUPER-KAMIOKANDE Collaboration], [arXiv:hep-ex/0209036]; D. Choudhury 

and A. Datta, arXiv:hep-ph/0606100. 
[5] F. Dydak et al, Phys. Lett. B 134, 281 (1984). 

[6] Y. Dcclais et al. [BUGEY Collaboration], Nucl. Phys. B 434, 503 (1995). 
[7] B. Armbruster et al. [KARMEN Collaboration], Phys. Rev. D 65, 112001 (2002) [arXiv:hep- 
ex/0203021]. 

[8] P. Astier et al. [NOMAD Collaboration], Phys. Lett. B 570, 19 (2003) [arXiv:hep-ex/0306037]. 
[9] A. Donini, M. Maltoni, D. Meloni, P. Migliozzi and F. Terranova, arXiv:0704.0388 [hep-ph]. 
[10] A. A. Aguilar-Arevalo et al. [The MiniBooNE Collaboration], Phys. Rev. Lett. 98, 231801 

(2007) [arXiv:0704.1500 [hep-ex]]. 
[11] A. Aguilar et al. [LSND Collaboration], Phys. Rev. D 64, 112007 (2001) [arXiv:hep- 
ex/0104049]. 

[12] M. Maltoni and T. Schwctz, arXiv:0705.0107 [hcp-ph]. 

[13] A. Kusenko, AIP Conf. Proc. 917, 58 (2007) [arXiv:hep-ph/0703116]. 

[14] G. C. McLaughlin, J. M. Fetter, A. B. Balantekin and G. M. Fuller, Phys. Rev. C 59, 
2873 (1999) [arXiv:astro-ph/9902106]; J. Fetter, G. C. McLaughlin, A. B. Balantekin and 
G. M. Fuller, Astropart. Phys. 18, 433 (2003) [arXiv:hep-ph/0205029]. 

26 



[15] J. Hidaka and G. M. Fuller, arXiv: 0706. 3886 [astro-ph]. 

[16] A. Kusenko, Int. J. Mod. Phys. D 13, 2065 (2004) [arXiv:astro-pli/0409521]. 

[17] T. Asaka and M. Shaposhnikov, Phys. Lett. B 620, 17 (2005) [arXiv:hep-ph/0505013]; 
T. Asaka, S. Blanchet and M. Shaposhnikov, Phys. Lett. B 631, 151 (2005) [arXiv:hep- 
ph/0503065]; N. Sahu and U. A. Yajnik, Phys. Lett. B 635, 11 (2006) [arXiv:hep-ph/0509285]. 

[18] S. Riemer-Sorensen, S. H. Hansen and K. Pedersen, Astrophys. J. 644, L33 (2006) 
[arXiv:astro-ph/0603661]. 

[19] F. Munyaneza and P. L. Biermann, arXiv:astro-ph/0609388. 

[20] S. Choubey, N. P. Harries and G. G. Ross, Phys. Rev. D 74, 053010 (2006) [arXiv:hep- 
ph/0605255]. 

[21] S. Choubey, N. P. Harries and G. G. Ross, arXiv:hep-ph/0703092. 
[22] R. L. Awasthi and S. Choubey, arXiv:0706.0399 [hep-ph]. 

[23] J. Hosaka et al. [Super-Kamiokande Collaboration], Phys. Rev. D 74, 032002 (2006) 
[arXiv:hep-ex/0604011]. 

[24] M. Apollonio et al. [CHOOZ Collaboration], Eur. Phys. J. C 27, 331 (2003) [arXiv:hep- 
ex/0301017]. 

[25] S. Goswami, A. Bandyopadhyay and S. Choubey, Nucl. Phys. Proc. Suppl. 143, 121 (2005) 

[arXiv:hcp-ph /0409224] . 

[26] E. K. Akhmedov, R. Johansson, M. Lindner, T. Ohlsson and T. Schwetz, JHEP 0404, 078 

(2004) [arXiv:hep-ph/0402175]. 
[27] A. M. Dziewonski and D. L. Anderson, Phys. Earth Planet. Interiors 25, 297 (1981). 
[28] C. H. Albright et al. [Neutrino Factory/Muon Collider Collaboration], arXiv:physics/0411123. 
[29] M. S. Athar et al. [INO Collaboration], INO-2006-01, 

http : / / www . imsc . res . in/~ino/OPENReports/INOReport . pdf . 
[30] P. Huber, M. Lindner and W. Winter, Nucl. Phys. B 645, 3 (2002) [arXiv:hep-ph/0204352]. 
[31] P. Huber, M. Lindner and W. Winter, Comput. Phys. Commun. 167, 195 (2005) [arXiv:hep- 

ph/0407333]. P. Huber, J. Kopp, M. Lindner, M. Rolinec and W. Winter, Comput. Phys. 

Commun. 177, 432 (2007) [arXiv:hep-ph/0701187]. 
[32] A. Ereditato and A. Rubbia, Nucl. Phys. Proc. Suppl. 154, 163 (2006) [arXiv:hep-ph/0509022]; 

A. Ereditato and A. Rubbia, Nucl. Phys. Proc. Suppl. 155, 233 (2006) [arXiv:hep-ph/0510131]. 
[33] S. Goswami and W. Rodejohann, arXiv:0706.1462 [hep-ph]; A. Bandyopadhyay and 



27 



S. Choubey, arXiv:0707.2481 [hep-ph]. 
[34] F. Ardellier et al. [Double Chooz Collaboration], arXiv:hep-ex/0606025. 
[35] X. Quo et al. [Daya Bay Collaboration], arXiv:liep-ex/0701029. 

[36] P. Huber, M. Lindner, M. Rolinec and W. Winter, Phys. Rev. D 74, 073003 (2006) [arXiv:hep- 
ph/0606119]. 



28 



